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Abstract. A number of scalar field models proposed as alternatives to the standard 
inflationary scenario involve contracting phases which precede the universe's present 
phase of expansion. An important question concerning such models is whether there 
are effects which could potentially distinguish them from purely expanding cosmologies. 
Vector perturbations have recently been considered in this context. At first order such 
perturbations are not supported by a scalar field. In this paper, therefore, we consider 
second order vector perturbations. We show that such perturbations are generated by 
first order scalar mode-mode couplings, and give an explicit expression for them. We 
compare the magnitude of vector perturbations produced in collapsing models with 
the corresponding amplitudes produced during inflation, using a number of suitable 
power-law solutions to model the inflationary and collapsing scenarios. We conclude 
that the ratios of the magnitudes of these perturbations depend on the details of the 
collapsing scenario as well as on how the hot big bang is recovered, but for certain 
cases could be large, growing with the duration of the collapse. 
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1. Introduction 

Recent years have witnessed tremendous advances in observational cosmology. High- 
precision data from observations of the Cosmic Microwave Background (CMB) and 
high redshift surveys have provided strong evidence for a nearly spatially-flat universe 
with a primordial spectrum of adiabatic, Gaussian and nearly scale-invariant density 
perturbations, in excellent agreement with the predictions of the simplest inflationary 
models [U E]. Despite these important successes, however, crucial questions remain. 
Paramount among these is whether the inflationary scenario can be embedded within a 
fundamental theory. While attempts in this direction continue, a number of alternative 
scenarios have been put forward motivated by string/M-theory. Among these are 
the pre-big bang [3j H] and the ekpyrotic/cyclic scenarios [HI El El E]- Despite 
their differences, an important ingredient shared by these models is the existence of a 
contracting epoch preceding a poorly understood bounce into the present expansionary 
phase of the universe. Interestingly, there are now many proposals for realising 
bouncing cosmologies within the context of string/M-theory [TOlCET], braneworld models 
P21 CL1 M, HH1 [33] and recently Loop Quantum Cosmology P0 US M, EU] , which add 
weight to the feasibility of scenarios of this type. Thus, given the successes of inflationary 
cosmology, on the one hand, and the importance of considering alternative scenarios, 
on the other, it is important to seek observational signatures that could distinguish 
between these alternative possibilities. 

An important tool in this connection is provided by cosmological perturbation 
theory. There are in general three types of perturbations, namely scalar, vector and 
tensor, which decouple at the linear level. A great deal of work has been done on 
studying the evolution of such perturbations in an expanding universe, both at the 
linear (for a review see [21]) and nonlinear levels (see e.g. [221 ESI [24J). A number of 
studies have also been made of the evolution of perturbations in collapsing models. In 
order to succeed as an alternative to inflation, any model of the early universe must 
provide an explanation of the primordial scalar perturbations, which are the seeds of 
structure in the universe. The amplitude and spectrum of these scalar perturbations are 
highly constrained by CMB observations, and although problems remain, both the pre- 
big bang and the cyclic /ekpyrotic scenarios have the potential to explain their origin. In 
this sense scalar perturbations may not be enough to distinguish between these scenarios. 

The next step is to consider other perturbation modes. Tensor perturbations 
are important as they can in principle be observed as a primordial gravitational 
wave background and can thus provide a way of distinguishing between different 
scenarios. Indeed, it has been shown that the cyclic /ekpyrotic scenario has very different 
predictions for tensor perturbations when compared to single field inflationary models, 
with the cyclic scenario producing a negligible level of gravitational waves with a blue 
tilted spectrum [25] . 

Until recently, however, little consideration was given to the evolution of vector 
perturbations. This is because these perturbations are generally assumed not to 
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be important in the inflationary scenario for two reasons. First, in an expanding 
FLRW universe first order metric vector perturbations decay and hence rapidly become 
insignificant. Secondly, once inflation has begun all matter except for the scalar field 
driving inflation is rapidly red shifted, with the result that the universe is effectively 
sourced solely by scalar field matter. This implies that vector perturbations play a 
minimal role since, as we shall see explicitly in the next section, a scalar field does not 
support vector perturbations at first order. 

In a collapsing scenario, however, the situation is potentially very different, as 
has been discussed recently in an interesting study [26J. Indeed, the first objection is 
no longer valid, since during a collapsing phase first order vector perturbations grow 
and in principal this growth could have important observational consequences for the 
collapsing scenarios [26]. The second objection is, however, still valid for the collapsing 
scenarios mentioned above, as these are also sourced solely by scalar field matter. Noting 
this, the authors of [26] proceeded by considering vector perturbations in presence of 
pressureless dust. However, the more natural setting to consider perturbations in this 
context remains that of scalar fields, where first order vector perturbations are absent. 
In view of this, it is important to ask whether nonlinear vector perturbations can be 
supported in such regimes, and how they evolve during collapsing phases. This is the 
aim of the present paper. 

We shall study the evolution of second order vector perturbations in a collapsing 
universe sourced by a scalar field. An important generic feature of nonlinear 
perturbations is that vector, tensor and scalar modes couple. Consequently, we expect 
second order vector perturbations to be produced, for example, by first order scalar 
mode-mode couplings. The analogous production of second order tensor perturbations 
has recently been studied [271 12H]- It is possible, therefore, that vector perturbations 
could provide a signature of a collapsing phase which would be absent or highly 
suppressed in an expanding universe. Indeed, the possibility that second order vector 
perturbations could act as seeds for large scale cosmic magnetic fields has previously 
been considered [29l [30] . as has their contribution to the polarisation of the CMB [31J. 

The structure of the paper is as follows. In Section [2] we introduce the flat FLRW 
model in the presence of a scalar field, and consider perturbations about this background; 
subsection 12.11 reviews first order scalar and vector perturbations, and an expression for 
the second order vector perturbations is derived in subsection 12.21 In Section [3] we 
consider a scalar field self-interacting through an exponential potential, which allows a 
power-law solution to the cosmological evolution. By considering different power-law 
regimes we are able to study concrete examples of interest analytically, and calculate 
and compare the amplitude of vector perturbations in each case in Section HI Finally, 
Section [5] contains our discussions and conclusions. 

Throughout lower case Latin indices take values 0, 1, 2, 3 and Greek indices 1, 2, 3. 
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2. Background model and Perturbations 

We consider as our background model a flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric in the form 

ds 2 = g ( ^dx a dx b = a{r) 2 {-dT 2 + 5 aP dx a dx p ) , (1) 

where r is the conformal time which is related to the coordinate time t through dt = adr. 
We take the universe to be sourced by a scalar field with a stress energy tensor given 
by 

Tab = <fi,a<fi,b ~ ^9ab(p' C (p,c ~ 9abV((f)) , (2) 

where V((f>) is the associated scalar field potential. The background evolution equations 
are given by the Friedmann equation 



8nG 



k2 



where H = a/a is the Hubble parameter and a dot denotes differentiation with respect 
to the coordinate time t, and the Klein-Gordon equation for the scalar field 

,5 + 3/^+^ = 0. (4) 

Throughout we shall use the formalism developed in [22] in order to give the vector 
perturbations up to second order in a flat FLRW universe sourced by a scalar field. 

In order to derive the perturbation equations we recall that the perturbed FLRW 
metric can, up to second order, be written in the usual form [22] 

g 00 = -a 2 (l + 2(A^ + A^)) , (5) 
9o* = -a 2 (BW + B^), (6) 

9ap = ^ (C + + O) > ( 7 ) 

where g^l is the background 3-metric and the superscripts (1) and (2) denote first and 
second order quantities, respectively. The perturbation variables can be decomposed, 
at each order i, as 

= /9g + , (8) 

C at3 = ( P KJ 9aP + l <a \p + G (a|/3) + C af3 I 9 ) 

with d a B ( a i] = d a C^ i] = d a C { ^ = C a ^ U) = 0. In this splitting, the variables <p® and 
ry{i) represent scalar perturbations, while vector and tensor perturbations are denoted 
by the superscripts (v) and (t), respectively. 

To proceed we require a gauge, which we choose to be the Poisson gauge, defined 
by /?(*) = T W = C ( a i] = 0. This gauge is a generalisation of the longitudinal gauge to 
include vector and tensor modes. 

We also note that we do not use the energy frame here. Instead, we follow [22] and 
use the normal frame, in which the scalar field stress-energy tensor can be identified 
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with the fluid stress-energy tensor in such a way that the velocity perturbations are 
related to the scalar field perturbations by 



t£> = - (10) 
a<p 

v^= -\ te ) +0«(0 (i) +0^ (i) )) . (ii) 

a<p v J 
We shall decompose the velocity perturbations as 

''aft 



with d a Va' ) = 0. We shall assume = throughout and, in the next section, we 

(vi) _ 



shall show that in our scalar field background we must also have v a = B\ =0. 
2.1. First order perturbations 

The first order perturbation equations in the present model have been considered 
extensively by many authors. An important feature of such perturbations is that the 
evolution equations for scalar, vector and tensor perturbations decouple and hence can 
be studied separately. 

Confining ourselves to the Poisson gauge, the metric perturbations and scalar 
perturbations (p^ 1 ' satisfy the relation = —ip^ 1 ' and their evolution equation becomes 

, m+ ^dy.rf +2 ^ H iy^ . (12) 

Employing the field equations together with the above fluid-scalar-field identification, 
to first order, the scalar field perturbation can also be given in terms of ip^ as 

47rG00 (1) = -(HtpW + (1) ) . (13) 

While considering first order scalar perturbations, we take the opportunity to define the 
curvature perturbations on comoving hypersurfaces 

3a 2 (l + w) ya'/a 3 , 

This quantity is extremely useful since for an expanding universe it is conserved on 
super-horizon scales. Thus it is the amplitude and spectrum of these perturbations 
which are important for the comparison of theoretical predictions with observations. 
We recall that the evolution equation for the first order vector perturbations is 
V 2^i) + 2HV 2 B^ = , 

which admits solutions proportional to 1/a 2 . 

Now, the right hand side of equation ( flOl) can be expressed as a gradient of a scalar, 
hence u« does not have a pure vector part, implying that a scalar field does not support 
a pure vector velocity perturbation at first order. We therefore have 

v™ = 0, (15) 
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in this case. Furthermore, the first order metric vector perturbations satisfy the 
momentum constraint equation 

= -8nG(p + p)v^ . (16) 

In line with common practice we shall employ Fourier decomposition of the 
perturbations. Considering Eq. ([TBI) together with Eq. (fl5j) implies that Ba 1 \k) = 
for k ^ 0, where Ba 1 \k) are the Fourier modes of the metric vector perturbations Ba ■ 
To proceed, therefore, we fix Ba = 0, so that at first order there are no metric vector 
perturbations present. 

Finally, to complete our first order analysis we comment on tensor perturbations. 
Since we are primarily interested in the evolution of second order vector perturbations, 
our interest in first order tensor perturbations is limited to how their couplings may 
produce vector perturbations at second order. However, their effects are likely to be 
subdominant as compared to the scalar terms, so we shall ignore the tensor perturbations 
at first order. Thus equations ([3])-(jI|) together with equation (Tl2l and the expression 
(fl3l) define a closed set of evolution equations for the background and the first order 
scalar perturbations in the Poisson gauge. 

2.2. Second order vector perturbations 

As was mentioned above, an important feature of perturbations at second and higher 
orders is that the evolution equations for scalar, vector and tensor modes couple. As 
a result, even in the absence of vector perturbations at first order, these modes can be 
generated at the second order by the scalar-scalar mode couplings as we shall see below. 

Following (22] (see also [23]), and using the Poisson gauge, the evolution equation 
for second order vector perturbations is given by 

_L V 2 (j g(, 2) + 2HBW) = MGVpUi + VpNl 



V^-VW^i + 8ttGTI?) , (17) 



where 



and 



1 I di) 1 



U aP = - 0gVS } - 7T<W«0 (1) ' 7 , (18) 



.(19) 



Taking V 2 of equation ( TT71) we obtain 

^V 2 V 2 (B^+2HB^) 

= AttG [V 2 0« VVff - 0? V 2 V 2 0« + ^v 2 ^ - <t>%v 2 ^~ 

+ 2vy;y^ - $ vyr + ^ . " ( 20 ) 
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This equation can be solved by giving </?W, which can be extracted from the first order 
evolution equation (fl2l) . together with the zeroth order equations. A non-trivial solution 
to ( TTT1) exists and has the structure 

C (x) 

B W = + (inhomo 

geneous terms) , (21) 

where C a (x) are arbitrary spatial functions and by "inhomogeneous terms" we mean 
the terms generated by the inhomogeneous part of Eq. ( 1201) . This solution has to be 
compatible with the momentum constraint equation, which in the Poisson gauge takes 
the particularly simple form: 

l V 2 i?^ = -8nGa(p + p){v^ - d a V- 2 d?vf ) , (22) 

where the second order velocity perturbations, v a , are given by Eq. ( TTTT) . We note that, 
in general, ^ and more importantly that Va 7^ 0. This is a crucial difference 
with respect to the first order case for which = forbids the existence of first order 
vector perturbations. 

Taking V 2 of the equation (1221) it can be rewritten using ffTTj) as 

j-V^B^ = (^g) + ^VV 1 ) - W 1 )^) - ffiW*) (23) 

- 47 rG(V 2 0«0« + - 0«V 2 0« - . 

The right hand side of this expression is in terms of first order scalar couplings which 
are in general non-zero. Furthermore, by direct substitution of Eq. (1231) into Eq. (120]) , 
we find that for (j2"U|) - (j2"T|) and (1231 to be made compatible we must have C a (x) = 0. 

This demonstrates that, in general, scalar-scalar mode couplings can produce vector 
perturbations at second order, even though they are absent at the first order. It also 
demonstrates that these perturbations are completely determined by the behaviour of 
the first order scalar perturbations, since the C(x)/a 2 part of the solution (1211) must be 
absent for a scalar field. 

It is also important to recall that there are a number of results suggesting that 
vorticity is zero in scalar field settings (see e.g. [221 123]). In our case, this is easy to 
demonstrate in the more familiar energy frame using the four velocity associated with 
the scalar field 

U a = J7<l>,a , 

where iV = \cf)£4>' b \ l l 2 ■, which implies 

w a t := h c [a h d b] u dtC = , 

where h a b = g a b + u a Ub- This is an exact statement and therefore valid at 
any perturbation order. So, we have a setting with non-zero second order vector 
perturbations and zero vorticity. 

Now, to determine the behaviour of second order vector perturbations concretely, 
we need to consider particular examples of contracting universes. We shall do this in 
the next section. 
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3. Analytic solutions 

In order to calculate the vector perturbations to second order we shall require solutions 
to the background and the first order equations. In general these can be obtained 
numerically. To proceed analytically, however, one needs to make assumptions. 
Concerning the background equations ©-(BJ, it is well know that there exists an exact 
solution if the field is self-interacting through an exponential potential [35J. The result 
is a power-law solution in which the scale factor grows with cosmic time as a = aot r . 
For these solutions H 2 oc <j) 2 oc V{4>) and hence the equation of state, w = - p = | 2 ~^|^ , 
is a constant. 

In both the standard inflationary scenario and scenarios involving a collapsing phase 
such as the ekpyrotic/ cyclic scenario, it is reasonable to assume that the equation of state 
will be approximately constant for a significant period of their evolution, and hence the 
power-law solutions are a very powerful tool for studying the dynamics of these scenarios 
analytically. Moreover, with the background evolving in accordance with a power-law 
solution, there is a well known analytic solution to Eq. (|T2l) . which then allows the first 
order scalar perturbations in this case to be determined analytically [36J. 

The form of the potential which gives rise to the power-law behaviour a = aot r is 

V(0)=y o e( _v ^)* ) (24) 

i 

assuming that </> is increasing. Changing to conformal time we have a = Oq~ v (— t) 13 ^, 



7i = -r—- and d>' = , 1 ^ t^--. In the following subsections we shall only be interested in 
cases for which r is negative and increasing towards zero, with — oo and zero representing 
the asymptotic past and future respectively. Using conformal time and taking Fourier 
transforms allows Eq. (fT2l) to be written as 

+ k 2 ^ l \k, r) +2(n'- H^j <P {1 \k, r) = , (25) 

where a prime denotes differentiation with respect to conformal time r and k = |k|. 
This equation can be solved in terms of Bessel functions. A common procedure is to use 
the transformation ip^\k, r) = (^) u(k, r) to rewrite the above equation in the form 

u{k, t)" + k 2 u(k, t) - — u(k, t) = , (26) 



x e , 

where 6 = TC/a<p', and hence ^- = r/(l — r) 2 r -2 . The solution to Eq. (1261) is then 
readily given by 

u(k, r) = (-r)2 (a k J u (-kr) + b k Y u (~kr)) , (27) 

where J u and Y v are Bessel functions of the first and second kind respectively, v — | | 
and a,k and b k are arbitrary constants. 

In all the scenarios we shall discuss, the perturbations represented by u(k, r) have 
their origin in quantum fluctuations which are normalised far inside the cosmological 
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horizon, but are pushed outside the horizon during the process of inflation or collapse 
[43J. Inside the horizon this corresponds to —kr — ► oo, and therefore in this limit we 
must have u(k, r) m 7^372 e~ lkr in order to match to the Minkowski vacuum. This limit 
allows the free constants to be fixed and we arrive at 



u(k,r) = ^e^ + ^V^H u (-kr) , (28) 
Ak 

where H v = J v (x) + iY u (x) is a Hankel function. The limit corresponding to the modes 
pushed outside the horizon is given by — kr — > 0. Expanding the Hankel function in 
this limit using 

*w=i^(i)"-f(^(ir +o ^- (29) 



7T 



we arrive at 



cp^(k,r) = C r k u -\-Ty-^-^ + D r k u+ \-r) u+ ^^ +0(r u+3 -^) 
+ B r k~ u+ \-r)- v+ i-^-r + A r k~ u ~\-T)- v - L i-^-r 

+ 0{t 3 -"-^), (30) 
where A r = fC r T(iy), B r = JC r jT(v> — 1), C r = fc>r V (J-i) > D r = — /C r ^"J^ 2 ) ' with 
/C r = -r^—K— { 2 u ~ie l ^ p+ ^ 1 5 —=A, — = ), are constants which are related to each other and 

depend on the initial data. This feature will play an important role in what follows. 

With the exact solution for ip^ 1 ' at hand, we are now able to evaluate the second 
order vector perturbations Ba 2 ^ from Eq. (|23j) . for a universe undergoing power-law 
evolution. A potentially important indicator to differentiate between collapsing models 
and expanding inflationary models could be sought in the ratio of the amplitudes of 
Ba 2 ^ between the expanding and contracting phases. 



4. Amplitudes of vector perturbations in expanding and contracting phases 

In this subsection we evaluate the amplitudes of Ba 2 ^ for the expanding and contracting 
phases. We restrict our attention to three scenarios which can be well modelled by the 
power-law solution introduced above. All these models can produce a nearly scale- 
invariant spectrum of first order scalar perturbations in a straightforward manner. The 
first order curvature perturbations ( must have a spectrum of this form in order to be 
compatible with observations. Thus the scale-invariance cannot be used to distinguish 
between these scenarios. We shall instead study whether the resulting amplitudes of 
the second order vector perturbations can be used to distinguish between these models. 

The first model we shall consider, r — > 00, corresponds to the standard expanding 
inflationary scenario. As we shall see, this limit of r gives rise to scale-invariant spectra 
for both <£>W and curvature perturbations. The second model we shall consider is 
the ekpyrotic/cyclic model, in which a contracting phase with r — > replaces the 
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inflationary epoch. This limit of r also gives rise to a scale-invariant spectrum for the 
(pW perturbations, but not for the curvature perturbations £. In the ekpyrotic/cyclic 
scenario part of the y^ 1 -* perturbation is matched onto the curvature perturbation when 
the universe undergoes a bounce, and hence the observationally important quantity ( 
has a scale- invariant spectrum after the bounce [371 EH]- Finally, we shall consider a 
dust-like contracting universe [391 BUI E] , which corresponds to a contracting universe 
with r — ► 2/3. While it has been pointed out that, unlike the previous two cases, this 
solution is not a dynamical attractor [12J S3], it is still of considerable interest since 
r = 2/3 gives rise to a scale-invariant spectrum for the curvature perturbations. We 
note that the limiting values of r we have discussed here are of course idealised. In a 
realistic setting r would approach but not reach such limits. This is important to keep 
in mind as the solutions for r = and r = oo are clearly ill defined. 

A comparison of these three cases is given in Ref. [43]. We shall now consider the 
three cases in turn, and evaluate Ba 2 ^ for each. We shall then proceed to calculate the 
ratios between the amplitudes of the second order perturbations Ba^ in the inflationary 
case and each of the contracting cases respectively. 

In what follows it is useful to define the constant quantities 

r 1 y/2r 

lr = i and Xr 



1-r ' " y/^Gl-r 
4-1. Expanding inflationary phase (r — > oo): 

This expanding phase is characterised by r — > oo (and v = \) [13]. Substituting in (J30D 
we obtain 

<p(k,r) = A^k-I + C^k" 1 * (-r) + (-r) 2 + 0(r 3 ). (31) 

For a perturbation to have a scale-invariant spectrum its Fourier components must be 
proportional to k~ 3 ^ 2 . Thus considering the expression above, we recover the well known 
result that as r — > 0, inflation produces a scale-invariant power spectrum for linear scalar 
perturbations. 

To evaluate the amplitude of the second order vector perturbations we consider 
the expression fl23|) . To proceed we shall employ the property that the Fourier 
transform (p(k,r) = f y?(x, r)e~ ik ' x <ix of the product <£i(x, T)(pi{x., r) is equal to the 
convolution ip^k^r) *^ 2 (^2,t) = / 0i(k 2 ,T)ip 2 (ki - fc 2 ,r)dk 2 - So e.g. f (V VvV ~ 
WVV)e~ Jk ' x dx = i j k|(ki - k 2 )(^A:2^fc 1 -fc 2 - 0k 2 l ^k 1 -k 2 )dk 2 . Then, by substituting 
( TSTT) into (1231) and using the Fourier transform we find 

7oo(7oo + 1) 



ktB^ = 4i ( 1 - ±J ) A^Coo 



too 



x J d 3 k 2 F(k 1 ,k 2 )(a k2 c kl ^ k2 - c k2 a kl _ k2 ) 

+ 0(r) , (32) 



where F(ki,k 2 ) = kf(ki - k 2 ) - [k 2 • (ki - k 2 )]k 2 , a k = k a ( b k = and c k = k~ 
with k — Ikl. 
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4-2. Recollapsing phase (r — > 0): 

This recollapsing phase is characterised by r — > (and v = |) [13]. Substituting in (1301) 
we obtain 

<^ (1) (fc,r) = Aok'l(-r)- 1 + C k~ l 2 + B ^(-t) + 0(r 2 ) . (33) 

Thus we recover the well known result that the linear scalar perturbations grow (diverge) 
in a collapsing phase if the so-called decaying modes are taken into account. This in 
turn implies (from (123]) ) that the amplitudes of second order vector perturbations also 
grow during a collapsing phase, as the inverse square of the conformal time. We can 
also see that ip^ has a scale-invariant spectrum in the r — > limit. 

In order to analyse this behaviour in more detail we substitute ( 1331) into (1231) and 
use the Fourier transform to obtain 

7o(7o - 1)\ AqC q 



kM v2) =a(i- 70(70 ~! ] ) 



r 2 



x J d 3 k 2 F(k 1 ,k 2 )(c k2 a kl _ k2 - a k2 c kl _ k2 ) 

+ 0(1). (34) 

It is interesting to note that the form of Eq. (123]) leads to a cancellation which results 
in the leading order time dependence in this case to be proportional to t~ 2 rather than 
r~ 3 as might be expected. 

4-3. Recollapsing phase (r = |j: 

This dust-like collapsing phase is characterised by r = 2/3 (and v = 5/2) [321 
Substituting in (!30l) we obtain 

<p {1 \k,r) = A2k-^(-Ty 5 + B2krl(-T)- 3 + 0(t--i) (35) 



Again we find that the scalar perturbations (p^ '(k,r) grow (diverge) in the collapsing 
approach to the singularity, but at a substantially different rate from the previous 
collapsing scenario. We shall see this has important consequences for the analysis of the 
next section. The second order vector perturbations in this case become 

4 7 ! - 17 7 2 + 75\ A2B 
jfci 4 BW = 8i 4 5 



AnGxl J (-r) 9 
3 ' 

d i k 2 F(k 1 ,k 2 )(b k2 a kl ^ k2 - a k2 b kl _ k2 ) 
+ 0(r~ 4 ) . (36) 

4-4- Comparison of expanding and recollapsing phases 

As discussed above, a potentially important signature of a collapsing phase can be 
obtained by comparing the ratio between the amplitudes of the second order metric 
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vector perturbations of the collapsing and expanding phases. In this subsection we shall 
evaluate this ratio for both contracting scenarios. 

Comparing the ekpyrotic collapsing (r — > 0) and the inflationary (r — > oo) cases we 
find the ratio of the amplitudes of the second order vector perturbations in these cases 
to be 

|B— | Ki A^ 1 +0(1)> 



| B («)«pj 'A^tI, 

where K 1 is a constant in time, which shows that this ratio grows with the duration of 
the collapse phase as l/r^ oll . 

Comparing the dust-like collapsing (r = 2/3) and the inflationary (r — > oo) cases 
we find the ratio of the amplitudes of the second order vector perturbations in these 
cases to be 

\T>(v2)coll\ A2B2 1 

K*1T7?-Z5- + 0(t£), (38) 



where K3 is a constant in time, which shows that the ratio in this case is radically 
different from the ekpyrotic case, growing with the duration of the collapse phase as 

Despite their usefulness these expressions are not, as they stand, sufficient to 
provide the complete information required for observational purposes. This is because 
we have not taken into account the constraints on the amplitude of the first order scalar 
perturbations. We shall proceed to implement these constraints in the next section. 



4-4- 1- Effects of imposing observational constraints from first order curvature 
perturbations: We have seen that in the scenarios which involve a collapsing phase the 
amplitude of the first order scalar perturbations (p^ grows as the collapse proceeds. At 
the end of the collapse, however, the first order perturbations must have the correct 
amplitude to be compatible with observations. It is therefore important to study 
the consequences of demanding that at the end of the collapse the first order scalar 
perturbations produced are equal to those obtained from observations. 

As was mentioned above, for expanding universes, and hence for inflation, the 
curvature perturbation on comoving hypersurfaces is conserved on super-horizon scales. 
This quantity is therefore used for comparison between theory and observations. In 
cosmologies with a collapsing phase we require the value of the curvature perturbation 
at the wavenumber corresponding to the largest scale on the CMB after the bounce to be 
equal to the required observational value, which is in turn equal to the value produced 
by a successful inflationary model. We note that since this study is only concerned with 
ratios we do not need to give the required observational value explicitly, and since the 
three scenarios considered here all produce the same spectral dependence we do not 
need to specify the wavenumber explicitly. 

For the scaling scenarios considered here we have 

a' _ ( r \ 1 a" _ ( r \ 1 

a \ 1 — r / r ' a' \1 — r J t 
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Thus in the expanding phase r^oowe find, using Eq. (1141) . that 
(exp = \— (C^k-i - B^k^-r) 2 + 0(r 3 )) , 



while in the collapsing case r^Owe obtain 

Ccoui = l—^ (—) (c AT* - 2B Q kh(-r) + 0{r 2 ) ) . (39) 



3 1 + w co m 

Finally in the r — > 2/3 collapse case we find 

Cco» 2 = (s.Ai-iC-r)- 8 + O(r-I)) . (40) 

3 1 + w coiZ2 V 3 / 

Considering the expression (j4"0l for ( co ii2, we recover the result that the collapsing 
scenario with a dust like equation of state produces scale-invariant curvature 
perturbations. 

In the ekpyrotic scenario (f^ is responsible for the curvature perturbation after the 
bounce. The argument used is that </?W an d ( mix at the bounce and part of tp^ at 
the end of the collapse is matched onto the curvature perturbations after the bounce. 
We therefore equate, at the lowest order, the perturbations 

(exp = Ipf colli , 

where tp represents the proportionality factor between the two phases in the ekpyrotic 
scenario [371 E] • This results in the constraint 

3 1 

Coo = r(l + w exp )A ip , 

£ Tcolll 

which can be used to obtain the ratio between the amplitudes of the second order vector 
perturbations in the collapsing and expanding phases from ( 1371) 

| B («2)coHl| / 4/3 



ib(-wi -M(irt^J +0{ cMl) - (41) 

This is a constant at the lowest order (independent of the length of the collapsing phase) 
thus indicating that the important factor determining the ratio is the proportionality 
factor ip in the ekpyrotic scenario, which is determined by the physics of the bounce. 

In the case of r — > 2/3 collapsing scenario, however, the curvature perturbations 
survive through the bounce and are thereafter conserved. Hence at the lowest order we 
can equate 

(exp (coll2 ) 

which implies 

1 + Wexp \ J_ 



0{r 2 cM2 ) . (42) 



Coo — B- , . 

3 \1 + W coU 2 J (-TT 

and from fl38|) 

| B («)coiI2| (l + U> col nV 1 



\ 1 + w exp J 



\Ti{v2)exp\ \ -I _|_ I 3 

i-> I \ 1 i- W PTT1 I l cM2 



Therefore, at the lowest order, the ratio of the amplitudes of the second order vector 
perturbations is proportional to t~J l , which increases as the collapse time tends to 0. 
Hence for a long collapsing phase we expect this ratio to become very large. 
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5. Conclusions 

A great deal of effort has recently gone into constructing alternative models to the 
standard inflationary scenario motivated by developments in string/M-theory. Many of 
these models involve a contracting phase. As a step towards distinguishing these models 
from the standard inflationary scenario we have studied the generation and evolution of 
vector perturbations in collapsing phases, in the context of scalar field cosmologies. 

There are a number of reasons why vector perturbations might provide useful 
signatures in models with contracting phases. Such perturbations are highly suppressed 
by inflation and are extremely difficult to generate at the early epochs after inflation. 
However, they would grow if produced during collapsing phases. 

Noting that first order vector perturbations cannot be supported by regimes purely 
sourced by a scalar field, we have considered second order vector perturbations. We 
have found that such perturbations can be generated by mode-mode couplings of the 
first order scalar perturbations and have derived their explicit dependence on first order 
perturbations. In principle our expressions allow the spectral dependence of second 
order vector perturbations to be determined, though in this work we have focused on 
their temporal behaviours in various early universe scenarios. Considering exponential 
potentials, which allow power-law solutions, we have studied the ratio of the amplitudes 
of second order vector perturbations in contracting and expanding phases. We have 
found that ignoring the details of how the hot big bang is recovered, the relative 
magnitudes of the second order vector perturbations in the two phases depend on the 
scaling solutions chosen. In particular we have considered two collapsing models, the first 
motivated by the ekpyrotic/cyclic scenario and the second by the dust-like collapsing 
scenario given in Ref. [39, 40J. 

For the first case, we have found the ratio to be independent of the length of collapse, 
at the lowest order, while depending on the details of the matching of the perturbations 
at the bounce. This is counter to the expectation that vector perturbations should 
be more prominent after a longer collapsing phase. This result is a consequence of 
the fact that the vector perturbations in this case grow as the square of the rate 
at which the scalar perturbations (p^ 1 ' grow, which in turn is a consequence of the 
cancellation mentioned in section 14. 21 In this case, cp^ 1 ' is responsible for the curvature 
perturbations after the bounce. Hence, once the observational constraint, namely that 
the curvature perturbations produced by this collapsing phase must have the same 
amplitude as the curvature perturbations produced by inflation, is imposed, the ratio 
of vector perturbations becomes fixed, independently of the collapse time. 

For the second case, we find that the ratio increases with the length of the collapse, 
being proportional to t co u~ 3 . This result is due to the fact that the second order vector 
perturbations grow more rapidly during the collapse than the square of the first order 
curvature perturbations. This implies that the magnitude of vector perturbation is no 
longer limited to being of the order of the square of the first order perturbations, but 
could grow to be larger. Hence, even when the constraint on the first order curvature 
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perturbations is imposed, there is still a dependence on the collapse time. This implies 
that the amount of vector perturbations present at the end of the collapsing phase could, 
in principal, be much larger than that present at the end of the inflationary epoch, since 
r — > as the collapse proceeds. 

It is important to note that the amplitude of the vector perturbations, Ba 2 \ 
calculated here is only valid if the power-law solutions of section [3] are applicable. 
These solutions are clearly only approximate and will not be valid through the entire 
inflationary or collapsing phases considered here. In particular, these behaviours must 
break down as the universe exits these phases and enters the hot big bang phase 
of evolution. In the inflationary scenario it is usual to assume that this occurs 
when the potential evolves such that it no longer gives rise to accelerated expansion, 
after which reheating occurs. In the collapsing scenarios the power-law evolution is 
expected to break down in the vicinity of the bounce, which is believed to lead to a 
radiation dominated phase. In principal, therefore, we would like to know how vector 
perturbations are affected by these transitions. This would, however, almost certainly 
require a numerical investigation, beyond the scope of the present study. 

In conclusion, a complete understanding of second order vector perturbations from 
a collapsing phase and their possible observational consequences would require a more 
detailed understanding of how such perturbations propagate between the scalar field 
phase and the hot big bang phase. However, by using power-law solutions we can 
compare vector perturbations produced in collapsing scenarios with those obtained in 
the inflationary scenario at a point immediately before the transition into the hot big 
bang phase. Our study of the two collapsing scenarios indicates that the observable 
differences between the collapsing models and the inflationary scenario could be large, 
particularly for the dust-like collapse, if we assume that the transition we have just 
discussed does not significantly alter the ratios we have calculated. 
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